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, 

A density-matrix method for the study of tagged states of neutral B 
Ph! mesons with arbitrary coherence properties is applied to several exam- 

r-| ! pies, including e + e~ production both at and above the T(45) resonance, 

and hadronic production. In the absence of coherence the only term mod- 
ulating the exponential decay of a neutral B meson behaves as cos Ami, 
while a sin Amt modulation is a signal of partial or full coherence. Decays 
to CP eigenstates are needed to fully specify the density matrix. We relate 
these results to more familiar expressions for the cases of the T(4S) and 
incoherent production. 



I Introduction 

Neutral B mesons undergo time-dependent oscillations with their antiparticles. This 
feature, first demonstrated for neutral kaons nearly half a century ago |T|], has been 
crucial in extracting fundamental information on the mechanism of CP violation from 
the decays of neutral B's. Moreover, the oscillations themselves have provided crucial 
information on the magnitude of electroweak couplings, and were one of the first pieces 
of evidence for a very heavy top quark 0. 

The oscillations are characterized by splittings Am between mass eigenstates. For 
the Bd = bd, the most recent world average is Am^ = 0.487 ± 0.014 ps -1 . For the 
B s = bs, only a lower limit 0, |J Am s > 15 ps _1 exists at present. 

1 To be submitted to Phys. Rev. D. 
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In the study of CP violation in decays of a neutral B meson, one frequently needs 
to know its flavor at the time of production. Was it a B° or a B ? Was it a B s or a 
B s l The dynamics of B meson production affords several methods for identifying this 
flavor. "Same-side" tagging methods []5], || [7|, |8| utilize the correlation of the flavor of 
a neutral B with the charge of a kaon or pion which is produced near it in phase space. 
"Opposite-side" methods utilize the associated production of bb in electromagnetic or 
strong interactions to tag a neutral B using the fragmentation products of the quark 
produced in association with it. The tagging methods are useful not only for the 
study of CP asymmetries, but also in the study of the oscillations themselves. For 
example, it is important to understand the systematic errors of tagging methods if a 
reliable estimate of Am s is to be achieved. 

The threshold for electromagnetic or strong production of a pair of nonstrange 
B mesons \M(B) = 5.28 GeV/c 2 ] is just below the T(4S) resonance, which lies at a 
center-of-mass energy of 10.58 GeV. At the T(AS), the reaction e + e~ — > BB produces 
the two mesons in an eigenstate of the charge conjugation operactor C, with eigenvalue 
rjc = —1. The flavor oscillations of neutral B's then manifest themselves as functions 
of the time difference t — i between their decays. Consequently, asymmetric e + e~ 
collisions have been adopted as a means of time-dilating the decays to enable the 
separation of their vertices 0, [TO ] . 

Another means of B production is through the decay Z — > bb, with subsequent 
fragmentation of the b or b to the neutral meson of interest. Here, the observed B 
and the tagging hadron (which could be any meson or baryon containing a b or b) are 
likely to be uncorrelated in their charge-conjugation properties, as has been assumed 



in several analyses (e.g., 11 



The question of coherence between the tagging hadron and the detected neutral B 
meson may not be so clear-cut in several cases intermediate between T(4S*) — > B°B 
(full coherence) and Z — > bb (little or no coherence). For example, just above the 
threshold for e + e~ — > B°B* or B B*°, if the photon in B* — ► B 7 or B*° — > B ^ is 
detected, the B°B pair will be in a state of r\ c = +1. If the photon is not detected, 
however, there may be an additional contribution from e + e _ — > B°B , in which the 
B°B pair has 77c = — 1. The relative probabilities P± of rjc = ±1 states are in any 
case unlikely to be equal. 

In hadronic bb production the subprocesses qq — > bb and gg — > bb generate a bb 
pair whose mass spectrum peaks at a scale of several times mj,. Additional bb pairs 
with an even sharper M(bb) peak near threshold arise from splitting of a virtual 
gluon: g* — ► bb. While incoherence has been assumed (e.g., ||12|| ) in analyses of flavor 
oscillations in hadronic B production, there may be effects of coherence if a B° and 
a B are produced in a state of low enough effective mass. This is particularly likely 
in the case of forward geometries (e.g., the HERA-B experiment at DESY |13| , the 
BTeV experiment at Fermilab H], and the LHCb experiment at CERN ||15|| ), in 
which the B and B are highly kinematically correlated. It is less likely to be the case 
for central B production (as in the CDF and DO experiments at Fermilab [T5, |T7| and 
the ATLAS and CMS experiments at CERN 



The practical effects of coherence cannot be ignored, since they lead to a char- 



2 



acteristic term proportional to e sin Amt in the dependence on proper time t of 
flavor oscillations. Such a term signals unequal probabilities for B°B production in 
eigenstates of positive and negative charge-conjugation eigenvalue. In the absence of 
coherence, the only terms present are proportional to e~ rt and e~ rt cos Amt. We thus 
advocate the inclusion of e _rt sinAmt terms in any analyses in which the presence 
of coherence is suspected. The present paper is devoted to the study of such effects. 
Although we shall generally speak of the (B°, B ) system, many of our results apply 
as well to neutral strange B mesons. 

In Section II we rederive a density-matrix formalism first introduced in J7|, |§ , using 
a more standard phase convention and correcting a sign error in the original references. 
This formalism is then applied to several cases, including mixed states with dilution 
of tagging efficiency (Section III), full coherence (Section IV), and intermediate cases 
(Section V). The means of fully specifying the density matrix is discussed in Section 
VI. The degree to which current and planned experiments can be expected to display 
coherence is given in Section VII, while Section VIII concludes. 



II Density-matrix description 



The density matrix is the appropriate means with which to discuss states with arbi- 
trary coherence properties. We work in a two-component "quasi-spin" space [2D, ^T| 
with initial basis states 



\B°) 



1 






1 



In this basis the most general density matrix p satisfying p 
written 
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-[1 + Q-a] 



(1) 



pt, Tr(p) = 1 can be 
(2) 



where Q describes polarization in quasispin space, Q 2 < 1, and <7j (i = 1, 2, 3) are 
the Pauli matrices. A pure state can be described by a density matrix with |Q| = 1, 
while a completely incoherent combination of B° and B with relative probabilities 
P B o and P-gO = 1 — P B o (a "mixed state" ) corresponds to a diagonal density matrix 
with Qi = Q 2 = 0, Q3 = 2P B — 1. One describes the density matrices for initial B° 
and B° by diag(l,0) and diag(0,l), respectively. 

The probability for a transition from an initial state denoted by the density matrix 
Pi to a final state denoted by p/ is then 



/(/) = Tr (p^pfT) 



(3) 



where T is the operator which time-evolves the state from i to /. Here / will denote 
an arbitrary coherent superposition of B° and B at time t, so that we shall be able 
to discuss decays to both flavor eigenstates (such as J/ipK*° — > J/tpK + 7r~) and CP 
eigenstates (such as J/i/)Ks,l)- The density matrix pf will take the appropriate form 
for each such final state. 
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It is most convenient to transform to the mass eigenstate basis 



\Bl) 



1 




\Bh) 




1 



(4) 



where "L" denotes "light" and "H" denotes "heavy," with the relation between mass 
and flavor eigenstates given by 



\B L ) = p\B°) + q\B") , \B H ) = p\B°) - q\B°) , \p\ 2 + \q\ 2 = 1 



(5) 



In a standard convention J22| one has q/p = e 2l/S , where (3 = Arg(—V* b V c d/V t lVtd), 
and Vij are elements of the Cabibbo-Kobayashi-Maskawa matrix specifying the charge- 



changing weak couplings of quarks. We then choose p = e 



1 3 



/\/2, 



q 



-0 



/V2- 



We shall neglect width differences in the following discussion. They are expected 
to be extremely small for nonstrange neutral S's, although they may be as large as 
10% for B s |23, [231]. In the mass eigenstate basis (|3J) the time evolution operator is 



(M'\T\M) = e 



-%M D t 



e- rt/2 diag(e- imL \ e ~ imHt ) 



-Ft /2 —imt ia 3 Amt /2 



(6) 



where m = (mu + m£)/2 and Am = tuh — tul- Transforming to the flavor basis ([J), 
we find 



(F'\T\F) = e - rt ^ 2 e- mt 



cos 

-2i0, 



Amt 



ie 2i(3 gin Amt 



Amt 
2 



COS 



Amt 



It is most convenient to express the density matrix in the mass basis as well: 
(M'\p\M) = ^2(M'\F'}(F'\p\F}(F\M} . 

F,F' 

We shall denote the density matrix in the mass basis by 

, 1 



P 



[1 + Q' ■ a] 



(7) 



(8) 



(9) 



the vector Q' in the mass basis is related to the vector Q in the flavor basis by 

Qi — Qs > 

Q> 2 = -(Q lS in2/3 + Q 2 cos2/3) , 
Q' 3 = Q 1 cos 2(3 - Q 2 sin 2f3 . 



(10) 

(11) 
(12) 



Since states which are pure B° or pure B correspond to Q s — ±1, Q\ = Q 2 = 0, 
their transformed density matrices are p\ = (1/2)(1 ± ax), respectively, since then 
Q[ = ±1, Q' 2 = Q'3 = 0. More generally, an incoherent state with Q3 = 2P B o — 1, 
Q 1 = Q 2 = corresponds to p' f = (1/2) [1 + (2P B o - 1)<ti]. 

The transition probability can now be written in terms of traces as 



1(f) = Tr (p^V /e 



-iM D t\ 



(13) 



For flavor eigenstates / corresponding to B° or B , let us take as examples J/ipK : 



*o 



J/%/jK + ti~ and J/ipK — > J/ipK~iT + , respectively. In the present convention both 
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Table I: Fractions f(H) of hadrons produced in b quark fragmentation. 



Hadron 


CDF (a) 


LEP 


Tf 


0.375 ± 0.023 


0.40 ±0.01 


B~ 


0.375 ±0.023 


0.40 ±0.01 


B s 


0.160 ±0.044 


0.097 ±0.012 


Baryons 


0.090 ±0.029 


0.104 ±0.017 



(a) Assuming equal fractions of charged and neutral nonstrange B mesons. 

decay amplitudes are equal to the same constant A, since they involve the quark 
subprocesses b — > ccs and b — > ccs, respectively. Then we find 

I ( oc e' rt [1 ± (Qi cos Amt + Q' 2 sin Amt)] . (14) 

The sign in front of the Q' 2 term was incorrectly stated in Refs. JF) and ||. 

As noted in Refs. @ and ||, the component Q' 3 does not appear in these expres- 
sions. We shall return to the question of its determination in Section VI. 



Ill Mixed state 

A mixed state of B° and 5° is one in which there are no amplitude correlations 
between the B° and B . Such a state will arise, in general, when a bb pair is produced 
with high enough effective mass that the b and b fragment independently. In this case 
we can consider a tagging method to indicate with probability P r the right-sign neutral 
B and with probability P w = 1 — P r the wrong-sign neutral B. The dilution factor 
V is V = P r - P w = 2P r - 1. 

Dilution can occur in various ways, depending on the tagging method. In opposite- 
side tagging at high M(bb), the opposite-side quark may fragment into a charged or 
neutral nonstrange B meson, a strange B meson, or a beauty baryon. These fractions 



have been measured at CDF [251 and LEP |31 and are summarized in Table I. 



The probability that a B° is detected as a B° is x^/[2(l + x 2 d )} ~ 0.18, where 



x d = Am d /T d = (0.487 ± 0.014 p S - 1 )(1.56 ± 0.04 ps) = 0.76 ± 0.03 @, g§. The 
corresponding probability of mis-detecting the flavor of a B s is very close to 1/2. 
Assuming that the other flavors are detected with unit probability, the CDF results 
imply P r ~ 0.85 and V = 0.70 while the LEP results imply P r ~ 0.88 and V = 0.76. 
In practice many other factors of course contribute to the dilution of a tagging method. 

Given a tag which should indicate the presence of a B° at time of production with 
tagging probability P B o, the corresponding density matrix elements in the mass basis 
are Q[ = 2P B o — 1, Q' 2 — Q' 3 — 0. The time- dependence of the flavor-specific final 
state / arising from either a B° or B decay is then given by 

( B° \ 

/ L « e~ rt [l ± (2P B o - 1) cos Amt] , (15) 
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without any sin Ami term. The quantity Pgo is usually determined empirically in a 
fit which also yields Am. 



IV Full coherence 

We now consider the case of fully coherent states of B° and B° produced in states of 
definite charge conjugation eigenvalue. We denote a C eigenstate of B u and B by 



=r0, ^ -f;0 



B u (p)B"(-p) + Vc B"(p)B»(-p) , (16) 



where r] C = ±1 is the eigenvalue of the C operator, and p and — p are unit vectors 
denoting the direction of the particles in their center-of-mass. The case of T(4S) — > 
B°B° corresponds to rjc — — 1. 

The states ^>c can be written in terms of mass eigenstates as 

*cfac = -l) = ^[B H (p)B L (-p) - B L (p)B H (-p)} , (17) 

^c(Vc = +l) = ^=[B L (p)B L (-p) - B H (p)B H {-p)\ . (18) 

These expressions allow us to write the elements Q\ of the density matrix in the mass- 
eigenstate representation and thereby to calculate the correlations between particles 
traveling along p (decaying at proper time t) and those traveling along — p (decaying 
at proper time t). We shall derive the results using both the one-particle formalism 
given in Section II and a two-particle formalism more suitable for joint distributions. 



A. One-particle description 

We consider for definiteness the case in which a flavor tag (B° or B°) is applied 
at a time t. Recalling that 

(B°\B L (-p)) = (B°\B H (-p)) = -L e * , (B°\B L (-p)) = -(B*\B H (-p)) = ±*r* 

(19) 

we can calculate the dependence on the tagging particle's decay time t of the pro- 
duction amplitudes of the states \B L>H (p)). For example, the state \B L {p)) ap- 
pears in (17) with coefficient —\Bh(—p))/V2, so it depends on t with coefficient 
-(\/2)e ifi e- Tl l 2 e- imHl for a B° tag and (l/2)e- if3 e - ri / 2 e- imHi for a B° tag. For states 
with both values of r\c we then find 

Vc = -1,B° tag 
f] C = -1,B° tag 



\Bl(p)) 
\Bh(p)) 

\B L (p)) 
\Bh(P)) 



1 if) -ri/2 
2 



2 



_ e -im H t 
e -irn L t 

e -irn H t ' 
e -im L t 



(20) 
(21) 
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Table II: Density matrix elements corresponding to correlated B°B° production in 
states of definite charge-conjugation eigenvalue rjc- 



Vc 


Tag (-p) 




Q' 2 


-1 


B° 


— cos Amt 


— sin Amt 




B° 


cos Ami 


sin Ami 


+1 


B° 


— cos Ami 


sin Ami 




B° 


cos Ami 


— sin Ami 



ric = +l,B°t&£ 



Vc = +1,5° tag: 



\Bl(P)) 
\Bh(p)) 

\Bl(P)} 
\Bh(P)) 



2 
2 



e ~im L t 
_ e -im H t 

e -im L t 
e -im H t 



(22) 



(23) 



Translating these pure states into normalized density matrices with unit trace, we 
find the results summarized in Table II. The component Q' 3 is zero. 

The density matrix elements in Table II can now be combined with the expression 
fll~4] ) to give joint rates for production of states with direction p decaying at time t 
and direction — p decaying at time i P7J. We find: 



I[B (t),B (i)]=I[B (t),B°(i)] 
I[B (t),B (i)]=I[B°(t),B (i)] 



- Tlt+F) [l-cosAm(t + T] C i)} 



(24) 
(25) 



The above expressions are consistent with those in the literature (e.g., [28]) and 
make physical sense. Their dependence on t + r/ci is mandated by Bose statistics. 
When t]c = — 1 and t = i, one never sees the decay products of neutral B mesons of 
the same flavor. This is also the case for t — i — when rjc = +1. When tags of each 
flavor are combined, the oscillatory terms cancel one another and one is left with a 
pure exponential ~ e _r< -* + ^. 

B. Two-particle description 

For entangled states such as described by in Eqs. (17) an d (18), one can use 
a direct-product notation [|29| . Our convention will be such that the first state in the 
direct product refers to the particle with direction p, while the second refers to that 
with —p. Typical direct products are then 



" 1 " 




' " 







1 






1 





" " 




' 1 " 




1 


1 

















(26) 



A spin-singlet state of two spin-1/2 particles is then represented by the four-component 



7 



vector 



IS = o> = -/=(!! -IT) 



1 

72 





-1 
1 





(27) 



and by the density matrix 



p(S = 0) = - 



In terms of the direct product representation, this can be written as 





" 








" 


1 





1 


-1 





2 





-1 


1 




















(28) 



p(S — 0) = -I <g> I - -<7i <g> <7i - -<r 2 <S> cr 2 - -cr 3 ® cr 3 
4 4 4 4 

where J is the 2x2 identity matrix. We have used the identities 



(29) 



01 <8> 0i 












1 " 










1 










1 










1 








_ 




" 1 










" 







-1 
















-1 
















1 



02 <X> 2 



-1 

1 

10 

-10 



I = 





1 






(30) 



03 ® 03 = 



Other useful identities are 



0i 09 02 = 



Writing the states and their time-evolution as in the previous subsection, we then 
find the corresponding density matrices 


1 





" 








—% 




" 








—% 










— % 













% 





(31) 





i 








, 02 <8> 0i = 





— i 








i 













% 














•c 



-1) = i e - r ('+*") 





-iAm(t-i) 



giAm(t— t) 







1 












I e -r(t+t) 



[I <g> I - a 3 <S> 3 



-(01 ® 01 + 2 
p(r/ c = +1) 



1 



r(t+t) 



_ e -iAm(t+F) 



E) 02 - 02 <8 
_ e iAm(t+t) 




1 



(32) 



01) sin Am(t — t)] 
= - A e~ T{t+T) [I ® I + 3 ® 03 



cr 2 ) cos Am(t — t) — (01 

1 







—(01 <S> 01 — 02 <S> 02) cos Am(t + t) + (o"i 

These results may be used to derive such expressions as (24) and (25) in an alternate 
way. 



02 + 02 <8> 0i ) sin Am(t + t)] 



(33) 
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V Intermediate cases 



The density matrices for r/c = +1 and i]c = — 1 can be added to one another. If the 
probability of states with rjc = ±1 is denoted by P± with P + + P = 1, the resulting 
elements for B°(p) production with a B°(—p) tag are 



B°(p), B°(-p) : Q[ = - cos Amt , Q' 2 = (P + - P_) sin Ami 



(34) 



with the signs of Q\ 2 changed for a B (—p) tag. The joint probabilities for production 
of (opposite,same) flavors of neutral B mesons at times t and t are then 



Opp 

Same 



j(t,f) = e~ r{t+F) [1 ± cos Amt cos Amt ± (P_ - P + ) sin Amt sin A 



mt] 



Equation (|35|) can be integrated with respect to time, with the result 



(35) 



dtl 



Opp 

Same 



1 ± 



ar 



cos Amt ± (P_ - P + 



1 + 



sin Amt 



(36) 

where a: = Am/T. As long as P_ ^ P + , the sin Amt term will be present. 

The two-particle description can also be applied to intermediate cases. For a state 
which is a mixture of rjc = +1 with probability P + and rjc = — 1 with probability 
P_ = 1 — P + , the density matrix is 



-r(t+t) 



P+ 




_p> g-iAm((+t) 





p_ 







.p) g-iAm(t-i) 
P_ 








p+ 



-r(t+t) 



7®/+ (P+ -P_)ct 3 ®ct 3 - 



e 



iAmt 



e 

-iAmt q 



iAmt 



HP- - p+ 



—e 



e 

-iAmt q 



Ami 



e 

-iAmt q 



iAmi 



(37) 



This last equation says, in particular, that in order to specify a one-particle state in 
which the matrix element Q' 3 is non-zero, one must not only have P + ^ P_, but the 
"tagging" particle (with decay time t) must also correspond to non-zero Q' 3 . 

These expressions hold for both nonstrange and strange neutral B mesons. They 
must be modified to take account of dilution effects such as those discussed in Section 
III. However, such effects should reduce the coefficients of the cos Amt and sin Amt 
terms by a common factor. 

Since x is very large for B s mesons, the presence of the sin Amt term may be 
difficult to demonstrate for them, unless one resolves the dependence on the "tagging" 
time t and does not integrate with respect to it. 
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VI Full specification of the density matrix 



As pointed out in Refs. [7[ || , it is necessary to observe decays to CP eigenstates and 
not just to flavor eigenstates in order to fully specify the density matrix, since the 
element Q 3 does not appear in any of the previous expressions for rates. We consider 
decays to J/ipK s and J/ipK L . 

Taking account of the negative CP of J/tpK s and positive CP of J/ipK Ll the 
decay amplitudes of interest are 



(J/ij;Ks\B ) = -(J/4>K S \B°) = (J/^K L \B°) = (J/^l|5°) = A' 
Then the density matrices for each final state in the flavor basis are 



Pj/^k s ^\A 



l\2 



PJ/ipK L 



/|2 



while in the mass-eigenstate basis they are 

1 



Pj/i>Kc 



Pj/i>K L 



A 



i\2 



A 



i\2 



1 - cos 2(3 
i sin 2(3 

1 + cos 2(3 
—i sin 2f3 



—i sin 2f3 
1 + cos 2(3 

i sin 2(3 
1 - cos 2(3 



(3* 



(39) 



(40) 



(41) 



We then recover the results of Refs. |7|, [8] , aside from a sign in the Q' 2 term which we 
correct here: 



I\J/1> 
l(j/^ 



K L 
K s 

K L 
K s 



oc e~ Tt { 1 ± [Q 3 cos 2(3 + {Q\ sin Amt - Q' 2 cos Ami) sin 2(3} 



oc e~ rt {1 ± [Q' 3 cos 2(3 - (Q\ sin Amt - Q' 2 cos Amt) sin 2(3} 



where / refers to a rate tagged with an opposite-side B, while / refers to a rate tagged 
with an opposite-side B. 

The determinations of Q' 3 and cos 2(3 are interrelated. Information on the sign of 
cos 2(3 would be useful in resolving the discrete ambiguity associated with extracting 
the value of (3 from that of sin 2(3 [|3D|. However, in eigenstates of C with t]c = ±1, the 
two-particle density matrix results indicate that the contributions from u 3 for the par- 
ticles decaying at times t and t are correlated. Thus, in order to prepare a state with 
Q' 3 7^ decaying at time t it appears that one must tag with a CP eigenstate decaying 
at time t. For example, in e+e" -> T(AS) -> B°B° -> (J/i(;K S;L )(p)(J/^K S;L )(-p), 
the effects of Q 3 7^ will always involve the term cos 2 2(3, so information on the sign 
of cos 2(3 is lost. 

An explicit calculation with the two-particle density matrix leads to the following 
time-dependent rates for a mixture of C eigenstates with probabilities P + and P_ = 

^-_[J/4,K S:L (t)J/^K SiL (i)] oc {1 + (P+ - P_) cos 2 2(3 
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- sin 2 2(3 [sin Amt sin Amt + (P_ — P + ) cos AmtAmi]} 
>[J/4K s , L (t)J/4K L , s (fi] ex {1 - (P+ - P_) cos 2 2/5 



(42) 



<it (it 

+ sin 2 2/5[sin Amt sin Amt + (P_ - P + ) cos AmtAmt]} . (43) 

As noted, (3 appears only through sin 2 2/5 and cos 2 2/5 = 1 — sin 2 2/3. 

In principle the reaction e + e~ — > T(45) — > P°P — ► (J/i(>Ks,l)(p)('K + 'X'~){— p) 
can provide additional information. The (normalized) density matrices for 7r + 7r~ 
production are 



in the flavor basis and 



Ptt+tt- 



1 

oil 



-17 



1 + cos 2a —i sin 2a 
i sin 2a 1 — cos 2a 



(44) 



(45) 



in the mass-eigenstate basis. (We have neglected the effect of penguin amplitudes 
here.) The use of Eq. (37) then leads to expressions involving the combinations 
cos 2(3 cos 2a and sin 2/5 sin 2a. 

In a state of definite charge-conjugation eigenvalue we find, by substituting the 
values in Table II and noting that Q' 3 = for such a state, that 



I [J/4 




oc e" r(i+i_) [1 =f sin 2/5 sin Am(t + r] C t)] 



I I J/4 ^ L ) oc e' v{t+l) [1 ± sin 2/5 sin Am(t + rj C t)} 



(46) 
(47) 



The results for rjc = — 1 agree with those in Ref. |28[, while for rjc = +1 the sign 
of the sin 2/5 sin Am(t + t) term is reversed. Again, although we have used the one- 
particle expressions based on the vector Q', these results can also be derived using 
the two-particle density matrices (32) and (33). 



VII Coherence expected in present and planned 
experiments 

The specific cases we have discussed so far range from fully coherent B°B° production 
at the T(4S") to uncorrelated production at high effective bb masses (as in Z° — > bb). 
A qualitative estimate of the degree of coherence expected between B° and B may 
be obtained by examining their relative orbital angular momenta. Since for a meson- 
antimeson pair with relative orbital angular momentum L the charge-conjugation 
eigenvalue is rj c = (— 1) L , the degree of coherence is expected to decrease as the 
accessible values of L increase. 

Suppose, for example, that the population of orbital angular momentum levels of 
B°B° is dictated by their statistical weights 2L + 1 up to a maximum L = L max . If 
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L max = the probability P + of rjc — 1 is 1, while the probability P_ of 77^ = — 1 is 0. 
If L max = 1 then P + = 1/4, P_ = 3/4. The general expressions are 

even ■ P ^ mffl ^~ ^ p L max (48) 

2(-^max +1) 2(L max + 1) 

r 11 r-> -^max L mSLX + 2 / .„-, 

L max odd : P + = — , P_ = — — . (49) 

^ V-^max + 1 21 max + 1 

Then one finds P + — P_ = (— l) Lmax / (L max + 1), and the magnitude of the coefficient 
of the sin Amt term in the time-dependence of a flavor eigenstate with a flavor tag 
decreases as 1/L max . 

A semiclassical argument can be used to estimate L max . Imagine a bb pair with 
squared cm. energy s to fragment into a pair of B mesons. The fragmentation process 
is limited to impact parameters b < 1 fm ~ 5 GeV -1 . Thus 

£max = kb , k = y/(s/4)-ml . (50) 

We now discuss the specific experimental cases mentioned in the Introduction, in 
decreasing order of likelihood of P°P° coherence. 

(1) Production at the T(AS) leads to a B°B pair in a state with L — 1, rjc — — 1. 
The full-coherence arguments of Section III apply. 

(2) Production through the process e + e~ — > P°P* + c.c. gives rise to a B°B pair 
with 77c- = +1, since B*° decays entirely to 7P . 

(3) Contamination of e + e~ — > B°B*° + c.c. fry e + e~ — > P°P° is likely. At the 
threshold for (vector + pseudoscalar) meson production, pseudoscalar meson pair 
production is also kinematically allowed, leading to some admixture of the 77c = — 1 
state. The relative fractions of the two processes are not well known but can be 
measured in principle. One can ensure against the i]c = —1 state by detecting the 46 
MeV photon in B*° decay. 

(4) Forward hadronic production of B°B will lead to a pair with effective mass 
typically no more that a few times the threshold mass of 2m&. The subprocesses 
qq — > bb and gg — > bb both favor low effective bb masses, while the gluon-splitting 
process g* — > bb will favor even lower effective masses. The corresponding value of 
k will then be of order m b , leading to L max = 0[m b ■ (1 fm)] ~ 25. Thus one might 
expect magnitudes of P_ — P + of at most a few percent. 

(5) Central hadronic production may lead to somewhat higher effective bb masses, 
especially if "opposite-side" tagging utilizes P's produced in the opposite hemisphere 
of the detector. The probability of B°B coherence is thus likely to be less than in 
forward geometries. 

(6) Production in Z — > bb is expected to lead to very little B°B° coherence, since 
k ~ M z /2 and L max consequently exceeds 200. 

The best prospect for studying the coherence effects we have mentioned here 
thus seems to be e + e~ collisions not far above the T(4S f ), where the rjc = +1 and 
i]c = —1 states are not necessarily equally populated. Ultimately, however, the 
question is an experimental one, and such effects can be studied at any energy and 
in any configuration by searching for the sin Amt term. 
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VIII Conclusions 



We have discussed the possibility of coherence of neutral B meson pairs, using a 
density-matrix approach which describes situations ranging from fully correlated pairs 
to mixed (uncorrelated) states. The density matrix is parametrized by a "polariza- 
tion" vector Q' describing a direction of "quasi-spin." Usual experiments determine 
only one component, Q[, of this vector, relating it to the dilution factor in flavor tag- 
ging. It gives rise to a characteristic modulation of exponential decay by a cos Amt 
term. In general there can appear a term proportional to sin Amt as well, which has 
not been taken into account in previous studies. This term arises from the compo- 
nent Q' 2 of the quasi-spin polarization vector, and is one signal of coherence. The 
component Q' 3 affects decays to CP eigenstates, and can be searched for by studying 
such final states as J/ipKs and J/i/)Kl- However, its investigation probably involves 
correlations between decays to pairs of CP eigenstates, and thus may require the 
production of a considerable number of B mesons. 

Acknowledgments 

We thank Yuval Grossman, John Jaros, Ady Mann, Yoram Rozen, and Sheldon 
Stone for helpful discussions. This work was supported in part by the United States 
Department of Energy through Grant No. DE FG02 90ER40560 and by the U. S. - 
Israel Binational Science Foundation through Grant No. 98-00237. J. L. R. wishes to 
thank the Technion - Israel Institute of Technology for gracious hospitality during 
part of this work. 

References 

[1] M. Gell-Mann and A. Pais, Phys. Rev. 97, 1387 (1955); A. Pais and O. Piccioni, 
Phys. Rev. 100, 1487 (1955). 

[2] J. Ellis, J. S. Hagelin, and S. Rudaz, Phys. Lett. B 192, 201 (1987); J. Ellis, J. 
S. Hagelin, S. Rudaz, and D. D. Wu, Nucl. Phys. B304, 205 (1988). 

[3] F. Palla, presented at Beauty-2000 (7th International Conference on B-Physics 
at Hadron Machines), Proceedings edited by P. Schlein, to be published in Nucl. 
Instr. Meth. 

[4] J. Jaros, presented at Beauty-2000 [||]. 

[5] A. Ali and F. Barreiro, Zeit. Phys. C 30, 635 (1986). 

[6] M. Gronau, A. Nippe, and J. L. Rosner, Phys. Rev. D 47, 1988 (1993). 

[7] M. Gronau and J. L. Rosner, Phys. Rev. Lett. 72, 195 (1994). 

[8] M. Gronau and J. L. Rosner, Phys. Rev. D 49, 254 (1994). 



13 



[9] D. Hitlin, presented at XXXth International Conference on High Energy Physics, 
Osaka, Japan, July 27 - August 2, 2000. 



[10] H. Aihara, presented at XXXth International Conference on High Energy 
Physics, Osaka, Japan, July 27 - August 2, 2000. 

[11] LEP and SLD Collaborations, D. Abbaneo et al, CERN report CERN-EP-2000- 
096 and SLAC report SLAC-PUB-8492, March, 2000 (unpublished), and refer- 
ences therein. 

[12] CDF Collaboration, F. Abe et al, Phys. Rev. D 59, 032001 (1999). 

[13] HERA-B Collaboration, presented by U. Uwer at Beauty-2000 @. 

[14] BTeV Collaboration, presented by T. Skwarnicki at Beauty-2000 |§. 

[15] LHCb Collaboration, presented by F. Muheim at Beauty-2000 0. 

[16] CDF Collaboration, presented by S. Bailey, J. Cranshaw, and M. Tanaka at 
Beauty-2000 [|. 

[17] DO Collaboration, presented by H. Evans and G. Gutierrez at Beauty-2000 [Q. 

[18] ATLAS Collaboration, presented by D. Rousseau at Beauty-2000 ||. 

[19] CMS Collaboration, presented by S. Banerjee at Beauty-2000 ||. 

[20] J. S. Bell and J. Steinberger, in Proceedings of the Oxford International Confer- 
ence on Elementary Particles, 19-25 September, 1965, edited by T. R. Walsh 
et al. (Rutherford High Energy Laboratory, Chilton, England, 1966), p. 195; 
T. D. Lee and C. S. Wu, Ann. Rev. Nucl. Part. Sci. 16, 511 (1966). 

[21] H. J. Lipkin, Phys. Lett. B 219, 474 (1989); Weizmann Institute preprint WIS- 
93-12, 1993 (unpublished). 

[22] I. Dunietz and J. L. Rosner, Phys. Rev. D 34, 1404 (1986); I. Dunietz, Ann. 
Phys. (N.Y.) 184, 350 (1988); I. Bigi and A. I. Sanda, Nucl. Phys. B281, 41 
(1987); I. Bigi, V. A. Khoze, N. G. Uraltsev, and A. I. Sanda, in CP Violation, 
edited by C. Jarlskog (World Scientific, Singapore, 1989), p. 175; Y. Nir and H. 
Quinn, Ann. Rev. Nucl. Part. Sci. 42, 211 (1992); J. L. Rosner, "CP Violation in 
B Decays," lectures given at TASI-2000, Boulder, CO, to be published by World 
Scientific. 

[23] D. Becirevic et al, preprint |hep-ph/0006135| (unpublished). 



[24] CPPACS Collaboration, S. Hashimoto et al, KEK preprint KEK-CP-101, |hep- 
| lat/0004022l (unpublished). 

[25] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 84, 1663 (2000). 



14 



[26] Particle Data Group, D. E. Groom et al, Eur. Phys. J. C 15, 1 (2000). 



[27] Here we have inserted an overall factor of e~ Tt to correspond to the exponential 
decay of the state with direction —p. It could just as well have been included in 
the (un-normalized) density matrix as a common factor. 

[28] I. Bigi et al. |22|. 

[29] See, e.g., Y. Ben-Aryeh, A. Mann, and B. C. Sanders, Found. Phys. 29, 1963 
(1999). 

[30] B. Kayser and L. Stodolsky, Max-Planck-Institut Report No. MPI-PHT-96-112, 
1996, |hep-ph/9610522l published in Beauty '97, edited by P. E. Schlein, Nucl. 
Instr. Meth. A , (1998); B. Kayser, National Science Foundation report NSF-PT- 
97-2, |hep-ph / 9709382j , in Proceedings of the Moriond Workshop on Electroweak 



Interactions and Unified Theories, Les Arcs, France, March, 1997; Y. Grossman 
and H. R. Quinn, Phys. Rev. D 56, 7259 (1997); H. R. Quinn, T. Schietinger, J. 
P. Silva, and A. E. Snyder, SLAC Report No. SLAC-PUB-8550, [hep-ph/0008021 
(unpublished). 



15 



